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Abstract 

Autoresonance is a phase locking phenomenon occurring in nonlin- 
ear oscillatory system, which is forced by oscillating perturbation. Many 
physical applications of the autoresonance are known in nonlinear physics. 

^U< The essence of the phenomenon is that the nonlinear oscillator selfadjusts 

to the varying external conditions so that it remains in resonance with 
the driver for a long time. This long time resonance leads to a strong 
increase in the response amplitude under weak driving perturbation. An 
analytic treatment of a simple mathematical model is done here by means 
of asymptotic analysis using a small driving parameter. The main result 

^^ is finding threshold for entering the autoresonance. 

in 

q ! 1 Introduction 

A hamiltonian system of two differential equations with a small parameter e 

v! — H v (u, v) = e/(r)cos(<^), v' + H u (u 7 v) = eg(r)cos(ip). (1) 



is considered. Here the right hand side represents a prescribed small external 
force which is fast oscillating: ip = $(r)/e. The /, <?, <&(r) are smooth functions 
of the slow time r = et. It is supposed that the unperturbed system (as e = 0) 
has the stable equilibrium point (0, 0) of center type in general position, which 
is taken as the initial point of the perturbed solution 
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u,«(t;e)|t=o = 0. (2) 

In this note we give an asymptotic solution of the problem (1),(2) as s — > 0, 
which is valid for large time t = 0(e^ 1 ). The main goal is to find conditions 

*The russian version of the paper is presented in "Matemat. Zametki" ("Mathematical 
Notes" in English) 



under which the system's energy H(u,v) + e[vf — ug] cos(^) grows up to the 
order of unity while the driver is being small: ef,sg = o(l), e— >0,. Such type 
phenomenon is referred to as autoresonance and one was studied by different 
authors [1-10]. 

Our construction is based on general two-parametric periodic solution of 
the unperturbed system: uo,v (t + to, E), (to £ R, E G (0, e ) C R) which 
is exists in a neighborhood of the equilibrium point (0,0). We set E — in 
the equilibrium solution wo,Uo(i, 0) = 0. The parameter E is a first integral 
of the unperturbed system, hence one can be identified as unperturbed energy 
E = H(u,v). 

We consider the general case, when both the period T = T(E) and the 
frequency to — 2tt/T(E) depend on the energy E so that T'(E) ^ 0. This 
dependence is a speciality of nonlinear system. It provides entering in autores- 
onance under different driver frequency. 

2 Anzatz 

An averaging method elaborated by Kuzmak [11] is here applied to asymptotic 
analysis of the solution. Sometimes this asymptotic approach is referred to as 
adiabatic approximation. The basic idea is that the leading order term of the 
asymptotic solution is taken as unperturbed solution 

:)-(^:§) + «(s&;l) + °c»- < t = e! » < 3 » 

under appropriate slow deformation of both the phase a = e~ 1 ^(t;£)/uj(E) 
and energy E = E(t,s). In the leading order term the problem is to find 
two functions ^,E(t;s) depending on slow time r = et. The main result is 
asymptotics of these functions. 

tt(T;e) = * (T)+e*i(T)+£ 2 * 2 (r) + O(£ 3 ), E(t,s) = £o(r)+e£i(r) + 0(e 2 ) 

(4) 
The initial energy is zero E(0, e) =0. If we find that E (r) = 0(1) as r = 0(1), 
it means that the energy of the system goes to order of unity for large times 

t = 0(e- 1 ). 

3 The slow deformation equation 

The property of T - periodic with respect to fast variable a identifies a class 
of functions u, v(a, r) which are used for asymptotic solution. This secular 
condition, as applied to the first order correction, leads to the equations for the 
desired slow varying functions. In the leading order term we obtain 

% = lo(E ), E' Q = ^J [d a u (a,E ),F( i p;T)}da (5) 



Here the integrand is skew-symmetric product of the c^uo = (u' (a, E),v' (a, E)) 
by the perturbation, which is F = (/, g) cos(ip) in the case of (1) . As usually in 
this way the relations (5) are looking as differential equations. But in the case 
under consideration we can convert theirs into algebraic equations. To this end 
one have to take into account that there is the fast variable (p — $(r)/£ with 
prescribed function 3>(t) in the original problem (1). It is not to hard to guess 
that in resonance case the desired fast variable a must to be related with the 
given if. Such type relation we refer to as 

Resonance requirement. Difference between driver frequency and free 
frequency is small 

*'-$' = O(e) (6) 

So we try to find asymptotic solution in the form (3), (4) under additional rela- 
tion (6). This approach may be realized only under some restrictions which are 
just autoresonance conditions. 

As the leading order term of the phase ^>q(t) is the given function $o( r ) 
hence the first equation in (5) is reading as an algebraic equation 

w(E ) = *'( T ) (7) 

for the energy Eq{t). Next, the second equation in (5) is reading as an algebraic 
equation for the first order correction of the phase function \Pi(r). In order 
to see it we have to take into account the resonance relation (6) as follows 
<p = uj{E)a - (* - $)/e = Lu{E)a - *i + 0{e). 

Analysis of higher order corrections (u n , v n ), n > 2 leads to linear algebraic 
equations for the corrections of both the energy and the phase function in the 
expansion (4). 

4 Results 

It is clear that in general case the nonlinear algebraic equations have no solu- 
tions. Requirements, under which they are solved, yield conditions of entering 
the autoresonance. 

Theorem 1 Let both uj'{E) ^ and the driver frequency coincide with the free 
frequency in the initial moment $'(0) = w(0). Then the equation (7) has unique 
growing solution Eq = Eq(t) for all r G [0, To] if and only if directions of the 
frequency variations coincide sgn<I>"(r) = sgau)'(E). 

Now the second equation in (5) is reduced to the algebraic equation for the 
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/ [u (T(E Q )s/2ir, E )), d s F(s - * i; r)] ds - $"(r)/u/(£ ) 
Jo 



(8) 



The equation can be solved only under specific conditions on the driver data 
F(s, t), $(t). In particular, one can see from (8) that the property: $"(0) = is 
need, because the unperturbed solution tends to zero: uq(T(Eo)s/2tt, Eq) — ► 



as energy tends to zero Eo(t) — > 0, r — ► 0. A crucial condition is imposed on 
the amplitude. In order to clear situation we consider a specific driver as given 
in (1). In that case the (8) is reduced to the trigonometric equation as follows 

$'(t)Kt) cos(tf i) - b(r) sin(*i)] = $"(t)/u/(£ (t)) (9) 

Here a, b arc Fourier coefficients depending on the unperturbed solution 

1 /■27T 1 /■27T 



... , [9(r)uo - f(r)vo] sin(s) ds, b= — [g{r)u - f(r)v ] cos(s) ds 
2tt Jo 27r Jo 

Thus the a, 6 are known as far as the unperturbed solution is known: (uo, vq)(o, E) 
under a = T(E)s/2ir, E — Eq(t). The following assertion is easily proved. 

Theorem 2 Let both the the right hand side in the original problem be given 
as in (1) and conditions of theorem 1 hold. Then equation (8) is solved if and 
only if the inequality 

a 2 (r) + b 2 {r) > [$"(t)/uu,'(E q (t))} 2 (10) 

holds. 

Condition (10) may be simplified. To this end an asymptotics of the unper- 
turbed solution as E — > can be used. In this way the following assertion is 
obtained: 

Corollar 2.1 Let both $"(0) = and the strong inequality 

^ 2 (0)[/ 2 (0)+ ff 2 (0)]>4|<D'"(0)M0)| 

hold in the initial moment. Then the equation (9) has a solution ty\ — ^i(r) 
on some finite interval r € [0, To]. 

5 Conclusion 

Requirements of the theorems 1,2 may be considered as conditions under which 
the system enters autoresonance, so that the energy of the system grows up to 
order of unity for large times t — r/e = C(e _1 ). In particular inequality (10) 
may be interpreted as threshold for entering the autoresonance. 

6 Restrictions 

For the first we have to remark that the formulas (3), (4) give only a formal 
asymptotic solution. The problem of justification of the asymptotics remains 
open. Moreover there remains an ambiguity in the the phase shift of the leading 
order term because of two roots of the equation (9). 

For the second we have to note that the formal asymptotic solution as given 
above is valid only for large times and one is not suitable on initial stage. 



In order to see this nuisance one have to analyze higher order terms in the 
expansion (4). They arc obtained from linear algebraic equations and in this 
way any additional requirements on the data do not arise. However higher order 
corrections have singularities as t — > 0, which became stronger as the number of 
the correction increases. For example, E\(t) = 0(t _1 ), ^(t) = C(r~ 3 )- Such 
type singularities indicate that the anzatz (3), (4) is not suitable on the initial 
stage, when t = et is small. Formulas (3), (4) give an asymptotic solution only 
for large times e~ 2 / 3 «K C(e _1 ). 

Asymptotic solution on the initial stage must be taken in another form in 
which the slow time 9 = e 2 / 3 i is used [12]. Such structure of the asymptotics is 
similar to boundary layer. Solution which is valid everywhere can be obtained 
by matching method [13]. 

7 Notes to bibliography 

Applications of the phenomenon similar to autoresonance was firstly suggested 
by Veksler [1,2], Andronov- Gorclik [3] and McMillan [4]. In particular, thresh- 
old for entering the autoresonance was pointed in [2]. Later, different systems 
were studied and entering the autoresonance was found in different cases, see for 
instance [5-10]. However accurate mathematical analysis of the autoresonance 
phenomenon does not undertake up to now. 
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